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ABSTRACT 

We report the discovery of a remarkable subpulse drift pattern in the relatively less studied 
wide profile pulsar, B0818— 41, using high sensitivity GMRT observations. We find simul- 
taneous occurrence of three drift regions with two different drift rates: an inner region with 
steeper apparent drift rate flanked on each side by a region of slower apparent drift rate. 
Furthermore, these closely spaced drift bands always maintain a constant phase relationship. 
Though these drift regions have significantly different values for the measured P2, the mea- 
sured P3 value is the same and equal to 18.3 Pi. We interpret the unique drift pattern of this 
pulsar as being created by the intersection of our line of sight (LOS) with two conal rings on 
the polar cap of a fairly aligned rotator (inclination angle a ~ 11°), with an "inner" LOS 
geometry (impact angle (3 ~ —5.4°). We argue that both the rings have the same values for 
the carousel rotation periodicity P4 and the number of sparks N sp . We find that N sp is 19- 
21 and show that it is very likely that, P4 is the same as the measured P3, making it a truly 
unique pulsar. We present results from simulations of the radiation pattern using the inferred 
parameters, that support our interpretations and reproduce the average profile as well as the 
observed features in the drift pattern quite well. 

Key words: Stars: neutron - stars: pulsars: general - stars: pulsar: individual: B0818— 41 



1 INTRODUCTION 



. ' The study of pulsars showing systematic subpulse drift patterns 
, provides important clues for the understanding of the unsolved 
5_l ■ problem of pulsar emission mechanism. Constraints provided by 
such observations can have far reaching implications for the theo- 
retical models, as exemplified by some of the recent results in this 
area (e.g. ( |Deshpande & Rankin (1999)| > and (Gupta et al. (2004)) ). 
In this context, pulsars with wide profiles — interpreted as, emis- 
sion coming from one magnetic pole of a highly aligned pulsar 
(pulsar with the magnetic axis almost parallel to the rotation axis) 
— can provide extra insights because of the presence of multi- 
ple drift bands, as illustrated in the recent studies of B0826— 34 
( (Gupta et al. (2004)) and (Esamdin et al. (2005)) ) and B0815+09 
( jOrao et al. (2004)] ?). 

B0818— 41 is a relatively less studied wide profile pul- 
sar with emission occurring for more than 180° of pulse lon- 
gitude. Discovered during the second Molonglo pulsar survey 
( (Manchester et al. (1978)) ), it has a period of 0.545 s and is rel- 
atively old, with a characteristic age of 4.57 x 10 8 years. The in- 
ferred dipolar magnetic field of this pulsar is 1.03 x 10 11 G, which 
is a typical value for slow pulsars. From a study of its average po- 
larization behaviour at 660 and 1440 MHz, (Qiao et al. (1995)) pre- 
dict that the pulsar must have a small inclination angle between the 
magnetic and rotation axes. 

We have carried out high sensitivity GMRT observations of 



B0818— 41, which bring out a unique pattern of subpulse drift, 
which is hitherto not reported. The main results from our initial 
observations and the interpretations thereof are reported in this pa- 
per. 



2 OBSERVATIONS AND DATA ANALYSIS 

We observed B08 1 8-41 at 325 MHz, on 24th February, 2003 at the 
GMRT, using the phased array mode ((Gupta et al. (2000))). The 
raw data were recorded with 0.512 msec sampling interval. During 
the offline analysis we further integrated the raw data to achieve 
the final time resolution of 2.048 msec. The duration of the obser- 
vations was about 31 minutes (i.e. 3414 pulses). 

The raw data were first dedispersed (with a DM value equal 
to 113.4 pc/cm 3 ) and then bad data points were filtered out from 
the dedispersed data. The final average pulse profile, illustrated in 
Fig. [T] was obtained by synchronous folding of this data with the 
Doppler corrected pulsar period. It shows a fairly steep sided and 
double peaked structure. The two components of the double peaked 
profile are separated by 86.6°±0.3°and are joined by a saddle re- 
gion. 

Single pulse sequences from our data reveal an interesting pat- 
tern of drifting. Fig.|2]shows a typical drift pattern over a sequence 
of 200 pulses. We see three different drift regions with two different 
drift rates — an inner region with steeper apparent drift rate flanked 
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Figure 1. Average profile of PSR B0818-41 at 325 MHz. The pulse am- 
plitude is in arbitrary units. 
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Figure 2. Gray scale plot of single pulse data (pulse # 200 to 400) of PSR 
B0818— 41 at 325 MHz, with the average profile shown on top. Signatures 
of radio frequency interference are present around pulse number 220, 298, 
338 and 397. 



on each side by a region of slower apparent drift rate. Multiple drift 
bands (typically 3 to 4) can be seen in the inner region for any sin- 
gle pulse, whereas the outer regions show only a single drift band. 
Though the regular drift patterns seen in Fig.[2]are quite common in 
our entire 3414 pulse sequence, there are regions of disturbed drift- 
ing, as well as significant sections of nulls. An illustrative 300 pulse 
sequence of this kind is shown in Fig. [3] Changes in drift rate, in- 
cluding curved drift bands (e.g. pulse # 720 to pulse # 780) are seen 
clearly. In this context this pulsar behaves similar to B0826— 34. In 
Fig. [3] some instances of the pulsar in null state are also seen — 
pulse sequences, 647 to 687, 693 to 705 and 780 to 835. 

We note that the inner and the outer drift regions are quite 
closely spaced with almost no discontinuity between the drift bands 
of the two region. Furthermore, the two drift regions are clearly 
locked to each other in phase — the subpulse emission from the 
inner drift region is in phase with that from the outer drift region 
on the right side, and at the same time the emission in inner drift 
region is out of phase with the outer drift region situated on the left 
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Figure 3. Gray scale plot of single pulse data (pulse # 600 to 900) of PSR 
B0818— 41 at 325 MHz, with the average profile shown on the top. Signa- 
tures of radio frequency interference are present around pulse number 625, 
720, 790 and 815. 
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Figure 4. The contour plot of the power spectrum of the flux as a function 
of pulse phase during a sequence of 200 pulses (pulse # 200—400). The left 
panel shows the power spectrum integrated over the entire pulse longitude. 
The upper panel shows the power integrated over fluctuation frequency. 



side. This phase locked relationship (hereafter PLR) is maintained 
for the entire stretch of the data and does not appear to get perturbed 
after intermittent nulling or during changes in the drift rate. 

We can characterise the observed drift pattern by estimating 
P™ (the measured time interval between the recurrence of succes- 
sive drift bands at a given pulse longitude) and Pf (the measured 
longitude separation between two adjacent drift bands). In practice, 
the P™ and P™ values need not correspond to the true values (P| 
and P 2 ') - see ( |Gupta et al. (2004)| >. 

To determine P™, we use the fluctuation spectrum analysis 
technique ((Backer (1970)1). Fig.|4]shows the phase resolved fluc- 
tuation spectrum for the 200 pulse sequence of Fig. [2] There is one 
strong peak at 18.3 ± 1.6 Pi (where Pi is the period of the pulsar). 
This feature is identified as P™ and the error bar in P™ is pro- 
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Figure 5. Same as Fig. 4, but for the inner drift region only 
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Figure 6. Auto correlation results for the inner region. Secondary maxima 
of the auto correlation function near the lag of 17.5 degrees. 



portional to the reciprocal of the total length of the pulse sequence. 
Fig. [5] shows the fluctuation spectrum for the inner drift region for 
the same pulse sequence. The P™ is found to be the same for both 
the inner and the outer drift regions, which is expected given the 
PLR. 

The P™ value for the inner drift region is calculated from the 
auto correlation function of the single pulses, averaged over the 
total number of pulses. The result is shown in Fig. [6] where the 
secondary peak is due to the correlation between the adjacent drift 
bands and the corresponding P™ value is 17.5°±1.3°. Because of 
the fact that one is not seeing simultaneous multiple subpulses from 
the outer drift region, P™ for this region can not be estimated using 
the correlation analysis method. Rough estimation of P™ value for 
the outer region — obtained from the slope of the drift bands in this 
region — is around 28°. Thus although the P™ values are the same, 
the P™ values are significantly different for the two drift regions. 



3 INTERPRETATION & MODELLING 

We interpret our observations and analysis results within the frame 
work of (Ruderman & Sutherland (1975)1 model and improved 
versions of the same (e.g. (Gil, Melikidze & Geppert (2003)| > and 
the references therein). According to these, subpulse drifting is pro- 
duced from a system of sub-beams (subpulse associated plasma 
columns). Sparks (sparking discharges within the vacuum gap), ro- 
tating around the magnetic axis under the action of an E x B drift, 
gives rise to a circulating pattern of sub-beams, and the time for 
one full circulation is referred to as the carousel rotation period, 
which we designate as P4. As pulsar radiation beams are widely 



believed to be arranged in concentric cones, it it natural to expect 
the circulating sparks to be distributed in annular rings on the po- 
lar cap (e.g. (Gil & Sendyk (2000)1), each ring giving rise to one 
cone in the nested cones of emission. To confirm and support our 
interpretations, we have carried out simulations of the expected ra- 
diation pattern for this pulsar. These are described at the end of this 
section, but their results are alluded to at the different stages of our 
interpretation, which are as follows. 

(1) The possibility that the observed drift pattern could be pro- 
duced from sparks circulating in one ring can be ruled out on the 
basis of the large difference in P™ values and the pattern of inten- 
sity distribution between the drift regions. A simple interpretation 
of the observed drift pattern is that the inner saddle region corre- 
sponds to an inner conal ring that is somewhat tangentially grazed 
by the observer's line of sight (LOS), while the outer drift regions 
are produced by the intersection of the LOS with an outer conal 
ring. The simulation results show that the inner drift region can be 
produced by the intersection of the LOS with three/four neighbour- 
ing sparks of the inner ring, whereas the outer drift region can be 
created by the intersection of the LOS with one spark of the outer 
ring, on each side of the inner ring. 

(2) The closely spaced inner and outer drift regions can be 
explained with a rapid transition of the LOS from one ring to the 
other. For the inner LOS geometry (negative /3), the LOS traverses 
the region between the inner and the outer region almost perpendic- 
ularly and as a consequence one can achieve quite closely spaced 
drift regions with reasonably well separated rings of emission. This 
is borne out by the simulation results (see Fig. [7}. 

(3) The observed PLR implies that the apparent angular drift 
rate is the same for both the inner and the outer rings. This can be 
achieved with two possibilities: different parts of polar cap plasma 
rotating at the same rate (i.e. the rotation is quasi-rigid), or differ- 
ent parts of the polar cap plasma circulating with different speeds 
which are fine tuned to maintain the apparent PLR. We do not see 
a natural way to achieve the later and hence, we suggest that the 
carousel rotation period, P4, and the number of sparks, N sp , are 
the same for both the inner and the outer ring. The simulations with 
two rings of emission with same P4 and N 3p , with 1 80° out of phase 
emission between the rings, successfully reproduces the PLR. 

(4) What are the likely values for N sp , P 3 (and hence P4), for 
this pulsar? To start with, we note that for a LOS grazing the inner 
ring, the P 2 m is expected to be close to P 2 *. So, N ap = 360/P 2 * ~ 
19 - 22. Now, for P 3 , there are two possibilities : (i) The P 3 m 
represents an unaliased true drift rate i.e., in successive pulses, the 
observer sees the same spark, shifted to a nearby longitude. In this 



case, Pi = P 3 m = 18.3 Pi and P 4 = N s . 



370 Pi 



(ii) The P™ is an aliased version of P 3 i.e., in successive pulses, 
a given spark drifts by an amount very close to the separation be- 
tween adjacent sparks (or a multiple thereof, depending on the or- 
der of the aliasing). As shown by (Gupta et al. (2004)), under the 
effect of aliasing, P 3 and P™ are related as, 



Pi 



(1) 



where k = INT[(n + l)/2] and / = mod(n,2) depend on the 
alias order, n. Note that we observe a negative drift rate for this 
pulsar. For an inner LOS geometry, this can happen for the case of 
an unaliased drift (as in case (i) above) — the ExB drift of the spark 
plasma is slower than the rotation of the star ((Ruder man (1976)) ). 
For an aliased drift, this can happen when n = 2, 4, 6 . . . (i.e k = 
1,2,3 ... ;l = 0). Taking the case for the lowest possible alias 
order (n = 2, k = 1, 1 = 0), and using P™ = 18.3 Pi, we obtain 
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Figure 7. Simulation of the subpulse drift pattern with simple dipolar ge- 
ometry for the case: a = -11°, P = -5.4°; drift rate= 20° /Pi. 



Pi = 18.3/19.3 Pi = 0.95 Pi and P 4 = N sp P 3 * = (18 - 
21) Pi m P™, This leads to the interesting result that the carousel 
rotation rate is the same as the measured P™. For this scenario, the 
entire spark pattern rotates with an angular velocity of 360° /P4 sa 
20° /Pi. 

To confirm these interpretations, we simulate the ra- 
diation from B0818— 41 using the method described in 
(Gil & Sendyk (2000) I. The goal of the simulations is to re- 
produce the following observed properties : (i) width of the pulse 
profile, (ii) the relative intensities of the drift regions and (iii) the 
overall pattern of drifting. Based on our interpretation, we use 19 
equi-spaced sparks in each of 2 concentric rings, rotating with a 
drift rate of 20° /Pi. The other free parameters in the simulations 
are the angle between the magnetic and rotation axes (a), the angle 
between the magnetic axis and the LOS (/3), the radii of the rings, 
the sizes and relative intensities of the sparks. We note that the 
presence of closely spaced drift regions with widely different P™ 
values, and our inference that the LOS intersects only one spark 
from the outer ring, all favour an inner LOS geometry (negative (3) 
over an outer LOS geometry (positive /3). 

For an inner geometry, we find satisfactory results for a fairly 
narrow range of values of a, (5: roughly from a — 12.2°and 
(3 = -6°to a = 7.5° and f5 = -3.35°. Fig. shows a typical 
simulation result, for the case of a = ll°and j3 = —5.4°. Here, 
the radius of the outer ring is 0.8 r p and the inner ring is 0.6 r p (r p 
being the polar cap radius). As can be seen, most of the observed 
features of the average profile and the drift bands are well repro- 
duced, including the significant change of P™ between the inner 
and the outer drift regions, which is clearly a geometrical effect. 
Interestingly, a 180° phase shift between the spark locations on the 
inner and outer rings is needed to match the PLR of the observed 
drift pattern (see Fig. [2}. For an outer LOS geometry, typical best 
fit results are for a = 2.8° and j3 = 7° . In this case, we find that for 
any reasonable geometry, (i) to reproduce the closely spaced inner 



and outer drift regions, the radial separation between the two rings 
needs to be rather small (0.72 and 0.82 r p ); (ii) to reproduce the ob- 
served P™ in the inner region (with the observed profile width), the 
sizes of the sparks need to be very small; (iii) furthermore, the wide 
difference in P™ values of the two regions is difficult to reproduce 
accurately. 

Finally, to investigate the relationship between P4 and P™, 
simulations have been tried using sparks of unequal intensities to 
bring out the presence of P4 explicitly in the fluctuation spectrum. 
Different values of P4 are achieved by using slightly different val- 
ues of the drift rate. In the general case, two different features can 
be seen in the spectrum, corresponding to P|" (caused by peri- 
odic drifting features) and P 4 (caused by the amplitude modulation 
due to the unequal sparks). For the appropriate choice of drift rate 
(equal to 20°/Pi), these two features merge/overlap in the spectrum 
to produce the measured P™ signal. For the real data, we believe 
that if the P4 is present in the spectrum, it is co-located with the 
aliased P3 signal, at P™ . 



4 DISCUSSIONS AND SUMMARY 

The remarkable subpulse drift pattern of B0818— 41 can be inter- 
preted as being due to an unaliased drift with a P4 ~ 370 Pi ~ 
200s, or an aliased drift with P 4 w 18.3 Pi « 10s. We 
believe that the latter is a more likely scenario for this pul- 
sar. An unaliased drift can not explain the occurrence of lon- 
gitude stationary sub-pulses and sense reversal of the drift rate 
that is seen in our data. For an aliased drift rate, this can be 
explained by small variations of the drift rate that move the 
P3 across the nearest Nyquist boundary ((Gupta et al. (2004) I). 
Though it is rare for the signature of carousel rotation period 
to be directly present as a low frequency feature in the fluc- 
tuation spectrum, we note that this has been claimed for at 
least two other cases : B0834+ 06 (([Asgekar & Deshpa hde (2005)) ) 
and B0943+10 in Q mode ((Asgekar & Deshpande (2001)f ~and 
l |Rankin & Suleymanova (2006)) ). 

It is interesting to compare our results with the the- 
oretical models for E x B drift, such as that given by 
( |Gil & Sendyk (20 00)). According to their equations (12) and (13), 
the complexity parameter a (which can be represented as Int(a) ~ 
2v + 1, with v being the putative number of cones), can be used to 
estimate the number of sparks in a ring as N sp — na = P4/P3. 
Taking P 4 = m Pi, we get : P3/P1 = m ir (2v + 1). Further, if 
P 4 = P 3 m , then from Eqn. 1 we obtain, P 3 */Pi = m (k m ± 1). 
Hence, n (2v +1) = (k m ± 1), which relates the number of cones 
to the measured periodicity and the alias order. Here, the plus sign 
applies for even alias orders, i.e. n — 2, 4, 6 • • •, I = and the 
minus sign for the odd alias orders, i.e. n = 1, 3, 5 • ■ •, I = 1. 
Thus, for our measured periodicity of 18.3 ± 1.6 Pi, the right hand 
side evaluates to 19.3 ± 1.6 for the inner LOS solution and first 
allowed alias order (n = 2). For the outer LOS solution, it eval- 
uates to 17.3 ± 1.6 for the first allowed alias order (n = 1). The 
left hand side evaluates to 15.7 ± 1.6 for v — 2 and to 21.9 ±1.6 
for v — 3. Thus, the outer LOS solution is more compatible with 
the pulsar having two cones, whereas the inner LOS solution is in 
closer agreement with a three cone model. It is quite likely that the 
LOS misses the third innermost cone in this pulsar. Note that drift 
with higher alias orders (e.g. n = 3,4) is incompatible with the 
above picture : either the number of sparks has to be twice as many 
(which is unlikely), or the P4 has to be shorter by a factor of 2. 

Furthermore, the complexity parameter can be used to es- 
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timate the screening factor r\ for the partially screened gap 
model, which describe the inner acceleration region in pul- 
sars ((Gil, Melikidze & Geppert (2003)) ). From equation (2) of 
( [Gil et al. (2006)) , P 4 /Pi" = a/2rj. For unaliased drifting 
(P 4 w 370 Pi), r] w 0.016 and 77 « 0.023 for v = 
2 and v — 3, respectively. For aliased drifting with P 4 = 
18.3 Pi, the corresponding estimates for 77 are 0.14 and 0.19. In 
(Gil, Melikidze & Geppert (2003)), calculated 77 values for differ- 
ent pulsars range from 0.032 to 0.36, with a typical value of ~ 0.2 
meaning that on the average the actual potential drop is of the order 
of 10% of the pure vacuum gap (( [Ruderman & Sutherland (1975)) . 
The 77 value for the aliased drift case for B0818— 14 is compara- 
ble to this, whereas for the unaliased drift case it will be the lowest 
amongst the known values of 77. This lends further support to our 
preference for aliased drift in the case of this pulsar. 

The attempt to determine n and P4 has been successful only 
for a few other pulsars : B0834+06 (P 4 = 14.8 Pi w 18.9 s, 
( |Asgekar & Deshpande (2005))), B0943+10 (P 4 = 37 Pi w 40 s, 
w=0, (Deshpande & Rankin (1999))), B08 09+74 (P 4 = 165 Pi w 
200 s,n = 0, (van Leeuwen et al. (2003)) ) and B0826-34 (P 4 w 
14 Pi w 25.9 s, n = 2 (Gupta et al. (2004)) ). We note that com- 
pared to these values, our results for B0818— 41 give P 4 ~ 200 s 
for the unaliased case and P 4 ~ 10 s for the aliased case. 

The permanent PLR between the inner and the outer drift re- 
gions is a unique feature in B0818— 41. It suggests that the rel- 
ative location of the circulating spark pattern in the two rings is 
strongly correlated, rather than being independent. Furthermore, 
the circulation rate is the same for both the rings, indicating a 
common electodynamic control in the entire polar cap. Arrange- 
ment of sparks located 180° phase shifted in the inner and the outer 
rings would correspond to the sparks that are maximally packed 
on the polar cap. Spark discharges occur in every place where the 
potential drop is high enough to ignite and develop pair produc- 
tion avalanche ( {Ruderman & Sutherland (1975)1), so sparks pop- 
ulate the polar cap as densely as possible. At the same time such 
arrangement of sparks should be quasi rigid, a property which is 
also suggested by the data. These results could be of significant im- 
plication for the physics of the polar cap. Multiple drift regions, 
with their drifts phase related to each other, are known in only 



a few pulsars, e.g. B0815+09 by (McLau ghlin et al. (2004)) , and 
B 1839-04 by jWeltevrede et al. (2006)) . The mirrored drift bands 
for B0815+09 are explained in the sparking gap model as emis- 
sion coming from an inner core gap (we refer as vacuum gap) and 
inner annular gap ((Qiao et al. (2004)1). The closely spaced inner 
and the outer drift regions can constrain the radial distances be- 
tween the rings of emission in the polar cap of B0818— 41. The 
radius of the inner and the outer ring from the simulation (Fig. |7J 



ent frequencies, are expected to provide new and interesting results. 
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pulsars by ( Gupta & Gangadhar a. (2003)) ) with a separate investi- 
gation method. 

To summarize, we find that B0818— 41 is a wide profile 
pulsar with unique drifting properties: closely spaced inner 
and outer drift regions, with apparently different drift rates, 
but permanently locked in phase to each other. To explain the 
observations, we find that the pulsar needs to be a fairly aligned 
rotator, with emission seen from two conal rings on the polar 
cap, likely viewed with an inner LOS geometry. Each ring has 
about 19 - 21 sparks, and it is very likely that the carousel rotation 
period is about 18.3 Pi ~ 10 s for this pulsar. The observations 
and analysis presented here constitute only a preview of the 
remarkable properties of B0818— 41. Follow up studies, including 
polarisation information and single pulse observations at differ- 
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